Mott Transitions in the 2D Half-Filled Hubbard Model: Correlator Projection 
Method with Projective Dynamical Mean-Field Approximation 
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The 2D half-filled Hubbard model is studied by a nonperturbative analytic theory of correla- 
tor projection. The dynamical mean-field approximation (DMFA) is reproduced at the first-order 
projection and then improved by systematic inclusion of spatial correlations at higher orders. A 
geometrical frustration induces a first-order Mott transition surface with a finite-temperature criti- 
cal end curve and related crossovers. Growth of antiferromagnetic correlations gives single-particle 
spectra strongly modified from DMFA with shadow bands and flat dispersions observed in high-Tc 
cuprates. 

PACS numbers: 71.30.-Hh, 71.10.Fd, 71.27.-Fa 



Analytic description of metal-insulator transitions 
(MIT) driven by electron correlation postulated by Mott^ 
has been one of the most intriguing issues in condensed- 
matter physics. This is also required for understand- 
ing high-Tc superconductivity^ that emerges around the 
MIT. In spite of substantial progress in the long history of 
studies'^"^, complete theoretical descriptions of electronic 
states around the MIT remain open. 

The dynamical mean- field (DMF) theory"^ unifies the 
Mott-Hubbard^ and Brinkman-Rice^ pictures within the 
Fermi liquid. In infinite dimensions where the theory 
becomes exact^, the first-order MIT line with a finite- 
temperature critical end point^^^ has been obtained for 
the Hubbard model on the fully frustrated Bethe lat- 
|-j(,gio,7_ jj-^ peaijg|;jQ dimensions, however, momentum de- 
pendence of the self-energy part and spin correlations 
ignored in the DMF approximation (DMFA) play other 
important roles near the MIT. Several studies have been 
done to restore such nonlocal correlations^^"^^. 

Much progress has been made on understandings of 
filling-control MIT in low dimensions by the operator 
projection method (OPM)^-'^. Based on the unique se- 
ries of Dyson equations or continued-fraction expan- 
sion^^ of the Green's function, the method systematically 
improves the mean-field, conscrving^^"^^ and Hubbard- 
type^'^*'^^ approximations. By taking account of spin 
correlations as well as the Hubbard-band splitting, this 
gives a unified scheme of the filling-control MIT beyond 
the early pictures^ This theory also reproduces nu- 
merically inferred key elements^''"^^ ; the diverging com- 
pressibility towards the MIT^", the four-band- like struc- 
ture^^ and flat dispersions around the momentum (tt, 0)^^ 
observed in high-Tc cuprates^"*. However, in particle-hole 
asymmetric cases, the self-consistent decoupling approx- 
imation^^ fails to reproduce a Mott insulator^^ due to 
insufficient estimates of local dynamics. 

In this letter, we propose an alternative method to 
Refs.^^"^'^ to take account of momentum dependences ig- 



nored in the DMFA. Combining the 0PM with a gen- 
eralized DMFA, we for the first time obtain an MIT 
phase diagram of the 2D half-filled Hubbard model in the 
parameter space of the local Coulomb repulsion U, the 
second-neighbor transfer t' and the temperature T scaled 
by the nearest-neighbor transfer t. t' introduces a geo- 
metrical frustration and a particle-hole asymmetry. Our 
results provide a comprehensive picture of the Mott tran- 
sition including the single-particle dynamics observed in 
insulating cuprates^^. 

At the first-order projection, our formalism hereafter 
called the correlator projection method (CPM) is reduced 
to the original DMF theory^. With higher-order projec- 
tions, spatial correlations are analytically restored in a 
systematic fashion. Here, the limit of the low energy 
and high momentum resolution is approached by simul- 
taneously including the lower-energy and longer-ranged 
spatial correlations along the same concept as the renor- 
malization group. This contrasts with the cluster meth- 
p^gi2,i3 -^jjgj-g ^j^g iiYmt of the high momentum resolution 
is taken after the low-energy limit. 

At the second-order projection, our projective DMFA 
gives the first-order MIT surface, critical end curve at 
T > 0, crossover between Mott- Hubbard and Slater in- 
sulators and metal-insulator crossover at higher T's in 
the {U, t' , T) space. Furthermore, beyond the original 
DMFA^^, antiferromagnetic (AF) or spin-singlet correla- 
tions produce the AF shadow structure and flat disper- 
sions observed in high-Tc cuprates^'''^^. 

This framework essentially satisfles the Fermi-liquid 
properties with a large Fermi surface and vanishing 
damping rates toward the Fermi energy and momenta for 
metal without any symmetry breaking, in contrast to pre- 
vious equation-of-motion approaches^'^*'^^. In contrast 
to cluster mcthods^^'^^, the CPM is free from cluster- 
size effects. It also gives an exact uniform solution of 
the soluble model with electron tranfer and exchange 
coupling between the nearest neighbors for SU{N) spins 
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with N — > oo^® in ID, while cluster methods give only 
approximate solutions^®. 

Wc consider the Hubbard Hamiltonian H = 

C/^^ 71x1^x4 on a square lattice, with the electron cre- 
ation, annihilation and number operators at a site x with 
a spin index s, c^.^, Cxs and rixs, respectively. 
and denote the summations over the nearest and 

next-nearest neighbors, respectively, is the chemical 
potential. In two dimensions, the magnetic phases do 
not appear at T > 0. The other symmetry breakings are 
assumed not to occur. 

The formalism of the CPM is based on the formal and 
hierarchical Dyson equations noted in the top panel [A] 
of Fig. 1. They are obtained by the projection in the 
equation of motion for Cxs^^- The explicit definitions for 
e^^'^^ and have already been given from the 
Hubbard parameters t, t' and U and cqual-timc corre- 

lations^^'^*. Here, e^^'^' simply yields the Hartree-Fock 
dispersions. 

To solve this continued-fraction expansion form re- 
quires an approximation to the highest-order self-energy 
part S„(w,k). Here, we restrict ourselves up to the 
second-order projection where n = 2. In earlier studies, 
S2 has been calculated by a two-site method'^" and a self- 
consistent decoupling approximation^^. However, with 
these methods, it is difBcult to discuss the MIT when a 
particle-hole asymmetry exists at half filling because of 
insufficient treatment of the local dynamics. Instead, wc 
employ a generalized DMFA to correctly calculate the lo- 
cal dynamics of S2(a;,k) by ignoring its k dependence. 
This is a systematic extension of the conventional DMFA 
since k) contains k dependences. 

The generalized DMF scheme is illustrated in the flow 
lines in the bottom panel [B] of Fig. 1: First, from an 
arbitrary S2,ioc, the local normalized self-energy part 
Gi,ioc(^) is calculated in the DMF procedure (1). Then, 
the normalized Weiss self-energy part Gi is obtained from 
Gijoc and £2,100- As in the procedure (2), the Weiss self- 
energy part Si directly obtained from Qi generates the 
Weiss Green's function Go by replacing Si with <Si in the 
Dyson equation. To calculate I]2,ioc(w), we employ the 
iterative perturbation scheme^ using Gai^^, k) in the pro- 
cedure (3). Then, the loop continues to (1) by replacing 
the original S2,ioc with S2,ioc obtained in the procedure 
(3) until the convergence of E2J0C is reached. A self- 
consistent solution is obtained after the iteration of this 
loop. The converged T,2{lo) is used to calculate G{uj,'k) 
by substituting into the [A]. We also note that in this 
loop, conservation laws for G{u),k) and 0(0;, k)Si(a), k) 
are satisfied. 

In the loop, Aecor/ {n){2 — {n)) is defined as the local 
part of £k This together with e^^^^) gives a source of 
the Hubbard band splitting, —tk and —ik are the non- 
local parts of e'^'^'' and s^^''^\ respectively, fk introduces 
a crucial k dependence of Si(a;,k) mainly through the 
superexchange interaction. The equal-time charge, spin 



and local-pair susceptibilities appearing in t^ may be de- 
termined self-consistcntly with the solutions of coupled 
correlator projection equations for two-particle operators. 
Here, however, we independently determine them from 
the two-particle self-consistent method^''. 

Now we show results from the above second-order 
CPM with projective DMFA. They reproduce impor- 
tant Mott-insulating character of single-particle spectra 
A{uj,'k) = — ;^Im G(cj,k), as shown in Fig. 2 for t' = 0: 
It produces a Mott gap A separating the two Hubbard 
bands where A seems to grow from with increasing U. 
At low T's, AF or singlet correlations yield AF shadow 
structure and flat dispersions around the (tt, 0) and (0, tt) 
momenta^^'^'^'^'*'^® in the Hubbard bands. Such signifi- 
cant modifications at the second-order CPM from the 
conventional DMFA^'' are due to inclusion of spatial cor- 
relations mainly through superexchange interaction. Be- 
sides, present results reproduce a direct gap for = at 
T/t = 0.02. k) shows a remarkable similarity to the 
QMC results^^. Deviations of the momentum distribu- 
tion at T/t = 0.02 from T = QMC data^^ turn out to 
be less than 10% for U/t = 4. 

Next, we discuss the MIT as a function of t' at half 
filling. The double occupancy (nfrii) exhibits a jump 
or a singularity depending on the first-order or second- 
order bandwidth-control MIT from arguments on the 
MIT in the effective action for spins and charges^'^*^. 
{n^ni) is plotted against \t' /t\ in Fig. 3. With decreas- 
ing T, the maximum slope of (n|n|) versus \t' /t\ in- 
creases. For U/t = 3 and 4, (n-^ni) exhibits a jump at 
\t'\ = t{^,[jrp([/, T) below the critical temperature Tcr{U). 
(Alternatively, if we vary U with t' being fixed, {n-^ni) 
exhibits a jump at some value UuiTit' ,T).) The jump 
indicates a first-order transition^. More correlated phase 
for \t'\ < t'yij^rj,{U,T) is insulating while less correlated 
phase with larger (n^ni) for \t'\ > tJ^^jrp(C7, T) is metal- 
lic. This is illustrated for the density of states p{uj) = 
if Ek ^(^' k) for U/t = 4 at T/t = 0.02 in Fig. 4. Here, 
an MIT occurs around t/ /t = —0.2525 with abrupt shifts 
in the spectral weights from the upper Hubbard band 
to the lower. The hysteresis seems to be suppressed; 
the higher-order projection enables us to escape from 
a metastable solution. Due to the particle-hole asym- 
metry introduced by t', the Kondo resonance does not 
appear in a particle-hole symmetric manner as in the 
fully frustrated Bethe lattice^. This first-order MIT be- 
low Tcr accompanies a discontinuous change in single- 
particle dispersions, while the dispersions continuously 
evolve at T > Tcr. The insulating dispersion with a 
choice oft = 0.25 eV agrees with that of La2Cu04^^. 

In this formalism, spatial correlations are restored as 
the k dependence of Si(w,k). Neglecting the k depen- 
dence of E„(ti),k) with finite n > 2 docs not automati- 
cally guarantee the Luttinger sum rule except the trivial 
case of n = 1: Including the k dependence of Ei(ci;, k) in 
the case of n > 2 modifies the shape of the Fermi surface 
and does not necessarily keep the Luttinger volume, if 
Si(a;,k) in the both limits of low energy and high mo- 
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mentum resolutions is not exactly obtained for finite n. 
At n = 1, the Fermi surface remains the same due to the 
absence of the k dependence of Si(a;,k). However, our 
metallic solutions satisfy the Luttinger sum rule with less 
than 16% deviations. 

As shown in Fig. 3, at low temperatures in the insu- 
lating region, (ni^ni) increases with decreasing tempera- 
ture for large U/t and small \t'/t\. Furthermore, {n^ni) 
decreases with increasing \t'/t\. Namely, (n-|-n|) seems 
to increase when the AF correlation is well developed^ 
This is interpreted from the following fact: With well de- 
veloped AF correlations, the coherence of quasiparticles 
(or in other words, coherent propagation of holons and 
doublons) increases, which results in the kinetic energy 
gain thus allows larger double occupancy. 

Fig. 5 shows the critical behaviors of the dou- 
ble occupancy; its jump at the first-order MIT 
A{nini) and the inverse of the maximum slope q = 
[|d(n-|-n|)/(i(f'/t)|]|(/|^f/ ([7,7-) obtained from Fig. 3. For 
fixed U, a critical point {t'^j.(U),Tcr{Uj) exists as a T > 
end point of the first-order MIT line (iMiTi^Mix)- The 
least square fits show that q diverges as (T — Tcr{U))~^-° 
for T > Tcr{U), while A{n^ni) vanishes as {Tcr{U) - 
y)0.23 fyj. J. ^ r„.(t/). Similar results have been ob- 
tained in infinite dimensions^''. 

Finally, the MIT phase diagram obtained for the 2D 
half-filled Hubbard model is shown in Fig. 6 with the 
first-order MIT surface (blue surface), critical end line 
(red curve), and metal- insulator crossover (green surface) 
that is determined from the temperature where the dou- 
ble occupancy varies most rapidly for fixed U and t'. Uct 
and Tcr vanish for t' = Q and increases with increasing \t'\ 
and/or T . Only at moderate t', the MIT phase diagram 
becomes similar to the DMFA results^' The first-order 
character seems to persist at T = due to absence of the 
particle-hole symmetry, unlike the particle-hole symmet- 
ric Bethe lattice"^. The MIT phase boundary at T ^ 
agrees with recent numerical results^^. 

In summary, a formalism to systematically include spa- 
tial correlations into the dynamical mean- field theory has 
been given. Mott transitions and the phase diagram of 
the 2D half-filled Hubbard model have been clarified by 
this correlator projection method. Local AF or singlet 
correlations prominently yield AF shadow bands, flat dis- 
persions around (tt, 0) and (0, tt) and a Mott insulator 
even at weak couplings, in agreement with the numeri- 
cal results^^"^^. In the phase diagram in the {U, \t'\,T) 
space, there exists a first-order MIT surface accompa- 
nied by a discontinuous change of single-particle disper- 
sions. The surface ends at a T > critical curve and 
a metal-insulator crossover appears at higher T's. The 
present theory qualitatively reproduces the dispersion in 
undoped cuprates. 

The work was supported by the Japan Society for 
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[A] Hierarchical 
Dyson equations 



G^(co,k): 



co-£'^'i)-Xi(co,k) 



-,(2,1) 



co-£(2'2)-22(«'k) 



[B] Generalized DMF scheme 

Second-order self-energy part ^2*^'^''*^) ~ 22,ioc('^) 



(1) Generalized DMF equation 

1 



Gl,loc(™)- atS^^ „(2,2) 

Normalized Weiss self-energy part 



(CO) 



(2) Weiss self-energy part = ^\ (to) 

Weiss Green's function ^(fS^M) - 



,1) 



(o-8'k'''-Si(«) 



(3) Iterative perturbation shceme 

22,loc(^) = 4Zrk,k.,q^o(^'k + q)^o(T,k'-q)^o(-T,k') 

k,k',q 

_ 4[fk'-^k'-q-fk+q+?'k+4ecor/(«)(2-(«))]^ 
^ k,k',q- («)(2-(«)) 



FIG. 2. (color) A density plot of the present results for 
A{ijj,\l) is shown in red color. Blue and black circles denote 
the previous results of the OPM with a self-consistent decou- 
pling approximation^^ and QMC results^'', respectively. The 
filled and the open symbols represent the momenta where the 
peak intensity is more and less than 10% of the maximum 
intensity, respectively. 




FIG. 3. (color) {n^ni) versus \t' /t\ at several U's and T's. 



FIG. 1. Schematic self-consistent loop of the currently pro- 
posed formalism of CPM. A'^ is the number of the sites. The 
CPM uses the generalized DMFA to calculate the local dy- 
namics of E„(tj,k), instead of the self-consistent decoupling 
approximation in the OPM^^. The first-order CPM is reduced 
to a conventional DMFA^. 
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FIG. 4. (color) Density of states p(tj) corresponding to 
metallic and Mott-insulating solutions. 
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FIG. 5. (color) The jump in {n-frii) (open symbols) and 
the inverse of its maximum slope q (filled symbols) ver- 
sus \t' /t\ are plotted against T/t with scaling functions. 
This gives a T > second-order critical point Tcr as an 
end point of the first-order MIT line {Tcr{U),t'^,{U)) with 

tuu) = t:{u,n,{u)). 




FIG. 6. (color) Phase diagram of the 2D half-filled Hub- 
bard model obtained by the present theory. The critical curve 
exists as illustrated by the red line. The green surface denoted 
as T, gives a metal-insulator crossover. 
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